Elastic Interactions of Cells 
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Biological cells in soft materials can be modeled as anisotropic force contraction dipoles. The 
corresponding elastic interaction potentials are long-ranged (~ l/r"' with distance r) and depend 
sensitively on elastic constants, geometry and cellular orientations. On elastic substrates, the elastic 
interaction is similar to that of electric quadrupoles in two dimensions and for dense systems leads 
to aggregation with herringbone order on a cellular scale. Free and clamped surfaces of samples of 
finite size introduce attractive and repulsive corrections, respectively, which vary on the macroscopic 
scale. Our theory predicts cell reorientation on stretched elastic substrates. 



Biological cells can exert strong physical forces on their 
surroundings. One example are fibroblasts, which are 
mechanically active cells found in connective tissue. In 
the early 1980s, Harris and coworkers found that fibrob- 
lasts exert much more force than needed for locomotion 
[Q. They suggested that strong fibroblast traction is 
needed in order to align the collagen fibers in the connec- 
tive tissue. Since cell locomotion is guided by collagen 
fibers, this results in a mechanical interaction of cells. 
The interplay of fiber alignment and cell locomotion has 
been analyzed theoretically in the framework of coupled 
transport equations for fiber and cell degrees of freedom 
[^. However, it is well known that cellular behavior is 
also affected by purely elastic effects, which were not con- 
sidered in these studies. For example, stationary cells 
plated on an elastic substrate which is cyclically stretched 
reorientate away from the stretching direction |^ , and lo- 
comoting cells on a strained elastic substrate reorientate 
in the strain direction [Q . Recent experiments show that 
adheringcells sense mechanical signals through focal ad- 
hesions p. In contrast to chemical diffusion fields, elastic 
effects are long-ranged and propagate quickly, and they 
are known to be important during development, wound 
healing, inflammation and metastasis j^. 

In this Letter, we consider theoretically the possibility 
of elastic interaction of cells. We focus on static forces, 
a situation which should apply to cells with restricted 
cytoskeletal regulation or to artificial cells which have 
a biomimetic contractile system without any regulation; 
the theoretical framework presented here for this case is 
a prerequisite for understanding the more complicated 
cases, e.g. the case of locomoting cells with a regulated 
response and dynamic force patterns In the static 
case, the elastic interaction of cells through their strain 
fields leads to forces and torques which can change their 
positions and orientations. If the cellular configuration 
can relax to equilibrium, the final configuration will be a 
minimum of the elastic energy. In the following, we derive 
the laws for elastic interactions of cells (which are mod- 
eled as anisotropic force contraction dipoles) and show 
how they depend on elastic constants, distance, cellular 
orientations, geometry and boundary conditions. 

If the distance between cells is much larger than their 



spatial extent, they can be modeled as point defects in 
an elastic medium. Elastic interactions of point defects 
have been discussed before for e.g. hydrogen in metal [Q, 
atoms adsorbed onto crystal surfaces and graphite 
intercalation compounds 1^. For each defect, the force is 
restricted to a small region of space, and the force distri- 
bution can be characterized by its force multipoles |IQ] 
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where f is the force density. The force monopole P is the 
overall force, which vanishes for inert particles. There- 
fore in the classical case, the first relevant term is the 
force dipole Pij, which describes the dilating/contracting 
action of the force distribution and has the dimension of 
an energy. Previous studies of elastic interactions of force 
multipoles were mostly concerned with isotropic force di- 
lation dipoles (that is Pij = P6ij with P > 0) and the 
finite sample size effect of free surfaces. The biological 
case which we discuss here is different in several respects. 
First, since cells can act as active walkers, there exists the 
possibility of force monopoles. Second, cellular force is 
based on actomyosin contractility and therefore leads to 
force contraction dipoles (that is P < 0). Third, adher- 
ing cells in most cases generate highly anisotropic force 
patterns, that is, the force dipole is not isotropic and 
will reorientate with respect to the surrounding strain 
field (e.g. on strained substrates). And fourth, in biolog- 
ical cases the elastic medium (e.g. the tissue) has clamped 
rather than free surfaces. In fact it is well known that 
cells become mechanically active only if their environ- 
ment can support enough stress, thus clamped boundary 
conditions are often needed to induce cellular activiation 

We assume that the elastic medium of interest (real 
or artificial tissues, elastic substrates) propagates strain 
like an isotropic elastic medium with a Young modulus 
in the order of kPa. For elastic substrates, the Pois- 
son ratio ly is close to 1/2 (incompressible case). For 
the following it is convenient to define A — X/ fi and 
c — 2)1 + X ~ /i(2 + A), where ^ and A are the Lame 
coefficients of the isotropic elastic medium. The incom- 
pressible case V = A/2{A+ 1) = 1/2 then corresponds to 
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the Hmit A ^ oo with A/c ^ Propagation of strain 
in an infinite isotropic elastic medium is described by the 
Green function jl^ 
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where r denotes the distance from the force center. For 
cells plated on an elastic substrate with a substrate thick- 
ness that is much larger than the displacements caused by 
cell traction, the relevant Green function is the one of an 
isotropic elastic halfspace with a free surface. Since such 
cells apply only tangential traction, we need to specify 
Gij only for the x-y- plane : 
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The elastic interaction energy W between two force dis- 
tributions can be written as a function of their force mul- 
tipoles and the Green function |]10| 
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where we sum over repeated indices and where indices 
after the comma represent partial derivatives. In the sec- 
ond line, the first line has been expanded twice and the 
definitions of Eq. have been used. 

The interaction between two force monopoles P and 
P' at r and r', respectively, follows from Eq. as 
W = —PiGij{r — J^')Pj- In the incompressible limit, this 
can be written as W = -(P ■ F' + (P ■ n)(P' • n))/87r/.tr, 
where n is the normalized separation vector between the 
two monopole locations (in three dimensions; in two di- 
mensions, an additional factor of 2 appears). This in- 
teraction is similar to the one between electric dipoles 
|p3| , thus we expect chaining to dominate large scale as- 
sembly, as confirmed by Monte Carlo simulations (not 
shown). However, since the force monopoles correspond 
to active movements, a model for cell locomotion is re- 
quired to fully treat this case. 

It is generally accepted that mechanically active cells 
exert only a very small overall force. Moreover, in most 
cases they are usually found to have highly polarized, 
that is, pinching force patterns iQ. In the following 
we therefore model cellular force patterns as anisotropic 
force contraction dipoles. The direction f of the pinch 
can be extracted from experimentally measured force 
patterns by determining the direction of the eigenvec- 
tor of the force dipole tensor corresponding to its largest 
eigenvalue. Then the force dipole tensor can be approx- 
imated as Pij — Pfifj. In many cases, the cell orien- 
tation following from the force pattern corresponds to 



the cell orientation following from overall cell shape or 
staining for actin fibers. For both locomoting and sta- 
tionary fibroblasts on elastic substrates, the magnitude 
of the force dipole can be estimated to be of the order 
of P « — 10~^^ J (this corresponds to a pinching pair of 
forces, separated by a distance of 60 ^m and each 200 
nN strong). The corresponding length scale (e.g. for dis- 
placements close to the cell) is (P/c)^/^ « 10 /im, which 
is somewhat smaller than a typical cell size (~ 50 /im). 
The interaction between two force dipoles Pu and P^^ at 
r and r', respectively, follows from Eq. (|4|) as 

iy(r,r') = -PHU.,z(i-,r') - PuG.jM^ - r')P^, (5) 

where u(r,r') is the displacement at r produced by the 
force dipole at r'. Since G ^ 1/r, the elastic interaction 
between force dipoles scales as ~ 

If the cells have isotropic force dipoles, their elastic 
interactions are well known: in infinite space, W — 
P^Gfj\j = [|lO| and an elastic interaction can only 
by induced by the boundary conditions H. On a semi- 
infinite space with free surface, W = P^Glf^^j = (2 + 
A)^P^/47r(l -I- A)cr'', thus the interaction is isotropic and 
repulsive Q. However, in most cases the cells will have 
highly anisotropic force dipoles. We start with the half- 
space and consider the following situation: one of the two 
interacting dipoles is fixed at the origin with vanishing 
polar angle. The other dipole is a distance r away with 
polar angle a. The polar angle of the separation vector 
is denoted by /3. Using Eq. (^ in Eq. (||), we find 

W{r. a, (3) = A_^J_/(a, p) (6) 
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with 



f{a, 13) =-[(4 + 3A) cos(2a) + 15A cos(2(a - 2/3)) 
8 

-l-(2-f A)(2-f 6cos(2(a-/3)) -f 6cos(2/3))] . (7) 

Depending on orientation, the interaction can be repul- 
sive or attractive. The attractive component leads to 
orientation dependent aggregation. 

In order to investigate this point in more detail, we 
consider force dipoles with a spherical hard core (cor- 
responding to a typical cell size). For A sa (vanishing 
Poisson ratio), the only favorable alignments will be side- 
by-side and the cells will assemble into linear strings, with 
their orientations perpendicular to the string direction. 
For larger A, aggregation will be much more compact. In 
the incompressible case, for a given angle f3 the optimal 
angle a follows as a,„j„ ~ 7r/2 -I- 2/3; the corresponding 
/ varies between -2 and -1.4. The optimal energy -2 is 
obtained for the four perpendicular configurations. Con- 
sidering only nearest neighbor interactions would lead to 
a square lattice at higher densities. However, at area 
densities beyond 7r/4, this structure is overpacked and a 
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herringbone structure results. Although the herringbone 
structure does not achieve the lower energy values of the 
square lattice, it can exist up to area density 7r/2-\/3 and 
will be favored for entropic reasons. For finite-sized clus- 
ters, surface reconstruction will take place. Moreover, 
for an increasing number of particles the interaction of 
Eq. leads to an increasingly rugged energy landscape 
with many local minimia due to the long range and ori- 
entation dependance of the potential. As a result of this 
metastability, unusual patterns like rings will form for 
certain initial conditions. In Fig. |l]we show some typical 
configurations. 

It is instructive to note that the oricntational part of 
the interaction in the incompressible limit is very simi- 
lar to that of linear electric quadrupoles in two dimen- 
sions (see Fig. p. This analogy is due to the fact that 
the corresponding interaction energy for linear electric 
quadrupoles, W = G'^jj^iPijPl^i (where G"='(r) - 1/r 
is the electric Green function), arises from contracting 
tensors of analogous symmetry. However, the near per- 
fect agreement in the angle-dependent part is an acciden- 
tal result for the incompressible case; moreover, electric 
quadrupoles interact with rather than with 

We now discuss the elastic interaction of anisotropic 
force contraction dipoles acting in a three-dimensional 
elastic medium, which follows by using Eq. (||) in Eq. (|^). 
This interaction is more complicated than the two- 
dimensional one, since its orientation dependance in- 
volves three rather than two different angles. A detailed 
discussion will be given elsewhere. Here we only discuss 
some high symmetry cases: for two parallel dipoles point- 
ing in z-direction and placed along the x-axis, we find 
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Thus this interaction changes sign as A varies through 
1 {v = 1/4). If one considers force dipoles arranged 
around a central dipole in the x-z-plane, similar consid- 
erations apply as in the two-dimensional case: for exam- 
ple, when considering only nearest neighbor interactions, 
a square arrangement of perpendicular dipoles with in- 
teraction energy W{r) = —(A + l)P^/47rcr^ is the most 
favorable one. However, if one now tries to continue this 
arrangement in the third dimension, frustration effects 
result. Depending on initial conditions, this enhances 
the occurance of irregular patterns. 

Due to the long-ranged nature of the elastic interac- 
tion, boundary effects will be very important 0,^. As 
an instructive example, we now discuss the elastic in- 
teraction of anisotropic force contraction dipoles in an 
isotropic elastic sphere of macroscopic radius R. For a 
free surface, this situation has been investigated before 
for both isotropic Q and anisotropic force dipoles [ p^ . 
However, no such treatment exists for clamped surfaces, 
which are expected to have larger biological relevance. 
For both free and clamped surfaces, one has to introduce 



image displacements, which can be determined using ex- 
pansions in vector spherical harmonics [ p^ . Again the 
general expressions will be given elsewhere, and here we 
only consider the high symmetry case of two dipoles both 
orientated in z-direction. Their direct interaction along 
the x-axis is given by Eq. (||). The image interaction fol- 
lows from inserting the image displacements of the first 
dipole (which are complicated functions expressed as vec- 
tor spherical harmonics) into Eq. (|^). For simplicity, here 
we report only the results for the imcomprcssiblc limit. 
We find 
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for a free surface and 
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for a clamped surface. Therefore free and clamped sur- 
faces introduce attractive and repulsive corrections, re- 
spectively. In Fig. H we show the interaction energies W 
for K = 2 [y — 1/3) as a function of distance x. We see 
that the image corrections can produce new minima in 
the full interaction potential. The main conclusion from 
Eqs. (§,|o|) is the fact that the image effects lead to cor- 
rections which operate on the macroscopic scale R. For 
the case of hydrogen in metal, this is known to lead to 
structure formation on a macroscopic scale [macroscopic 
modes) ||7|,p^. In the biological case, the boundary in- 
duced pattern formation competes with structure forma- 
tion on cellular and elastic scales, which results from the 
direct elastic interaction. Therefore we expect that in 
a theory for cellular densities modes in elastic media of 
finite size, hierarchical structures will result. Like for hy- 
drogen in metal, image effects can be also expected to 
lead to incoherent deformations (fracture). 

The simplest case of the interaction of a cellular force 
pattern with an elastic strain field is the case of a single 
cell plated on an elastic substrate which is homogeneously 
stretched along the x-direction by applying stress p > 
at the sides. Then {ux,Uy) = (j)/E)(x, —lyy), with E be- 
ing the Young modulus. The interaction energy follows as 
W — —PiiUij = —{Pp/ E){cos^ a — z^sin^a), where a is 
the polar angle describing cell orientation. Since P < 0, 
the cell will reorientate perpendicular to the direction of 
stretching {a = 90°). In the case of compression, p < 0, 
the cell will reorientate parallel to the direction of com- 
pression (a = 0°). Although the dynamic aspects of elas- 
tic interactions of cells are out of the scope of the present 
work, it is worth noting that on cyclically stretched elas- 
tic substrates, p periodically changes sign and in order to 
maintain a stationary state, the cell might try to avoid 
both tensile and compressive strain. One easily calcu- 
lates that this corresponds to a = arccos + v). 
For V K, 0.4, this yields a ~ 60°, in excellent agreement 
with experiments 
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FIG. 1. Typical configurations for anisotropic force dipoles on an incompressible elastic substrate, (a) Minimal energy 
configuration for two dipoles (f = -2). (b) Linear electric quadrupoles have the same optimal configuration, (c) - (f) Small 
clusters are subject to surface reconstruction. For five particles, (e) and (f) are nearly degenerate (f — -11.59 and f = -11.43, 
respectively), (g) At high densities, herringbone order results, (h) In terms of energy, the square lattice is most favorable. 
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FIG. 2. Elastic interaction energy W in units of /cR^ for two parallel anisotropic force dipoles of magnitude P in an 
isotropic elastic sphere with radius R and elastic constants c and A = 2 (Poisson ratio u — 1/3). For this value of A, the direct 
elastic interaction is repulsive (dotted lines), (a) Free surface: the image correction (dashed line) is attractive and generates a 
new minimum in the full interaction potential (solid line), (b) Clamped surface: the image correction is repulsive. 
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